Statistical description of eigenfunctions 
in chaotic and weakly disordered systems beyond universality 
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We present a semiclassical approach to eigenfunction statistics in chaotic and weakly disordered 
quantum systems which goes beyond Random Matrix Theory, supersymmetry techniques, and ex- 
isting semiclassical methods. The approach is based on a generalization of Berry's Random Wave 
Model, combined with a consistent semiclassical representation of spatial two-point correlations. 
We derive closed expressions for arbitrary wavefunction averages in terms of universal coefficients 
and sums over classical paths, which contain, besides the supersymmetry results, novel oscillatory 
contributions. Their physical relevance is demonstrated in the context of Coulomb blockade physics. 
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Since Chladni's famous experiments on vibrating 
plates two centuries ago, the morphology of eigenfunc- 
tions in wave (mechanical) systems has aroused curios- 
ity 0. In particular, the statistical description of eigen- 
functions in quantum (or wave) systems with diffusive 
or chaotic classical (or ray) dynamics has been an inten- 
sive research field for more than 20 years |2j ■ Besides the 
pure academic interest, the study of the spatial prop- 
erties of irregular eigenfunctions provides specific theo- 
retical input for experimentally relevant quantities in a 
wide range of disciplines 1] including optics, acoustics, 
oceanography |3|, quantum chaos, and mesoscopic con- 
densed matter In the latter, for instance, the energy 
levels of interacting electrons in quantum dots are influ- 
enced by wave function correlations entering into inter- 
action matrix elements H, IS EL 

In the semiclassical limit |8j, characterized by typi- 
cal classical actions S c \ much larger than Planck's con- 
stant, i.e. h e s = h/S c \ *C 1, statistical measures for the 
spatial structure of eigenfunctions can be deduced from 
purely classical quantities, though in a highly non-trivial 
way. The most spectacular example of such a connec- 
tion is embodied in Berry's Random Wave Model (RWM) 
eigenfunctions of classically chaotic systems possess 
the same statistical properties as Gaussian random fields 
with a universal spatial two-point correlation function. 

A theory beyond the universal results of the RWM has 
been rigorously derived for disordered systems by an ex- 
act mapping of the quantum problem onto a supersym- 
metric field theory, the Nonlinear Sigma Model [6( . Semi- 
classical in spirit, this Diffusive Sigma Model expresses all 
results in terms of the classical diffusion propagator, and 
its success has motivated efforts to extend such methods 
to clean chaotic systems. The results of such (still conjec- 
tured) Ballistic Sigma Models (BoM) [lOH 

are obtained 

by replacing the diffusion propagator for the disordered, 
metallic regime by a suitable ballistic counterpart. Two 
features of the BoM render calculations difficult, already 
at the level of the first non-universal contribution beyond 
the RWM (i) the statistical distribution of eigen- 



functions is not Gaussian, in clear contrast to the RWM, 
and (ii) different eigenfunctions are not independent. 

In this Letter we present an alternative approach to 
eigenfunction statistics in which different eigenfunctions 
are described as independent Gaussian fields [12j. Em- 
ploying simple Gaussian integrations, we provide closed 
expressions for general wave function averages for chaotic 
systems in the semiclassical regime, in terms of universal 
coefficients and sums over classical paths. We demon- 
strate that an arbitrary average contains, besides the 
universal (RWM) result, a system-dependent contribu- 
tion (e.g. from confinement) composed of a smooth (di- 
agonal) and an oscillatory part. We resolve the apparent 
contradiction to the supersymmetry techniques by show- 
ing how the results of the BoM beyond universality are 
retained from our result by neglecting all oscillatory con- 
tributions. This finding for chaotic dynamics is used to 
extend our method to the metallic regime of disordered 
systems, while the physical relevance of the oscillatory 
terms found is illustrated for Coulomb Blockade conduc- 
tance peak statistics. 

Defining statistical averages. Consider a set of d 
normalized, real solutions ip na (n a £ {ni, . . . , n,j}) of 
the Schrodinger equation with non-degenerate, generally 
non-consecutive eigenvalues E ni < . . . < E nd . We study 
expressions of the general form F(ip ni , . . . , ip nd ) where 
each entry is a vector ijj na — (i>n a (ri), ■ ■ ■ , VW^L )) 
with f a components depending on different positions 



rf,i = 1, 
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by convention). Upon varying n a inside each window 
[n a — N/2, fi a +N/2], while keeping all differences n a — np 
fixed, the function F(ip ni , . . . ,tpn d ) will exhibit fluctua- 
tions. The spectral average (indicated by caligraphic let- 
ters) of the function F is then naturally defined as 



N/2 



(1) 



--N/2 



It is essential to keep in mind that T is fluctuating itself, 



2 



depending on the size and location of the energy windows 
and on the set of positions rf . In the semiclassical limit 
one has n a ^S> N ^> 1 for all a. 

Averages playing a key role are the spatial two-point 
correlation matrices R Q with entries defined as 



f 



N/2 



R «=N 2^ ^ a+ s(r7)^ a+ s(rJ) . (2) 

s = -N/2 

Local Gaussian conjecture. A local Gaussian theory 
for the spectral averages is based on two assumptions: 

(i) ipn a -N/2i ■ ■ ■ > 4>n a +N/2 are N realizations of an f a - 
dimensional random vector, denoted by v a , with distribu- 
tion P a (v a ) = (27r) _/ °/ 2 \/detRQ :L exp [—^v a .R~ 1 v a ] ; 

(ii) v a and vp are independent random vectors for a ^ (3. 
In the local Gaussian theory all averages can be expressed 
through the correlation matrices R Q by means of 
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(vi, . . .,Vd)Pi(vi) ■ . .P d (vd)dvi . ..dv d ■ (3) 

The local Gaussian conjecture states that T = J- G , 
namely that any two eigenf unctions of a classically 
chaotic system behave like two independent Gaussian 
random fields, each of them uniquely characterized by the 
exact two-point correlation functions. This is a natural 
generalization of the ideas presented in 0, H3. fljl fl6| 
with far reaching consequences when a subsequent semi- 
classical approximation is consistently used. 

Semiclassical expansion. In the semiclassical regime 
each correlation matrix is expanded into a constant part, 
a leading-order (in 7i e ff) fluctuating part, and higher- 
order terms [l7j]: Ra = A- 1 l a +A~ 1 'R a -\0(fi^), where l a 
is a f a x f a unit matrix and A the system area for the 2d 
case considered. The matrix ^4 _1 I Q defines a Gaussian 
distribution P RMT (v a ) = (A/27r)^/ 2 exp [-{A/2)v a -v a ] 
yielding the Random Matrix Theory (RMT) results 
■pRMT £ or q,.^ state. It is straight forward to fac- 
torize each probability distribution as 



Pn 



:P RMT (v n 



exp 



• (Io+Rq)" 



det(I Q + R Q ) 



(4) 



As a key step we note that the semiclassical approach 
consistently keeps terms up to second order in the fluc- 
tuating part of the correlation matrices, since the semi- 
classical approximation neglects terms of order 0(fi.^ 2 ) 
@ while R a ~ 0{h\£). Taylor expansion of the proba- 
bility distributions in Eq. Q to second order in R Q and 
substitution into Eq. © yields 



RMT 
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with universal (system-independent) coefficients 



T 
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,, - {F{v 1 ,...,v d )q iij {v a )) RMT , 



Fi,j,k,i = (F(vi,...,v d )qij,k,i(v a )) 



RMT 
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(F(vi, . ■ . ,Vd)qi,j(v a )q kt i{v/3)} 



RMT 



Here (. . .) RMT denotes an average with respect to the 

distribution Yl v =i P RMT (^i)^ anc ^ 9$) are simple poly- 
nomials of the components of its argument v: 



H,j(v) = ^(Av lVj -Sij) , 
q t ,j,k,i(v) = ^[qi,j(v)qk,l{v)+2q i ,k(v)qj,l( i ")] ( 7 ) 
4 

Semiclassical correlation function. The fluctuating 
part of the correlation matrix is obtained from the semi- 
classical Green function and is given by [l4l 



Rti j = Si.i + Jo(fc„Jrf ~fj 

1/2 /T 

li,3 
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as a sum over classical trajectories joining rf with rf, 
at fixed energy E ria ~ n a A, where A denotes the (con- 
stant) mean level spacing. The Bessel function Jq(x) 
with kn a — y / 2mE ria /h is the contribution from the 
unique direct trajectory between rf and rf. The sum 
in Eq. JSJ is taken over all non-direct orbits jij with 
actions S 7M {rf , r^, E Ha ) — J p- df and semiclassical 
prefactors -D™. . where the stability and topology of each 
path enters |§|. The window function T(x) = sin(x)/x 
suppresses contributions from trajectories with traversal 
time T li . = dS li /dE na larger than the characteristic 
time tat = 2H/NA related to the energy average JTJ. 

Equation (JjJJ, supplemented by the definitions I|tjl7|) 
and Eq. (jHJ) for R 1 ^ , is our main result. We illustrate its 
power by computing the intensity distribution I(ij)(ri 
S(t— Aip(ri) 2 ), a prominent and frequently studied 
measure for wave function statistics. While the evalua- 
tion based on supersymmetry methods is quite involved 
using our Gaussian approach (with d = 1, /i = 1) we 
trivially find for the universal polynomials qi,i{v) = 
\{Av 2 - 1) and gi.i,!,!^) = \{Av 2 - l) 2 - \A 2 v A . The 
coefficients ^ are then obtained without any further in- 
tegration, and Eq. JSJ yields the intensity distribution 



l G (t) = 



pi,i 

l + ^-(*-l) + 



(3-m+t 2 ) 



in terms of closed orbits (through R 1 ' 1 ) starting and end- 
ing at r\. Equation Q includes both the universal limit 
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e~i /y/2irt (the Porter-Thomas distribution Q) and the 
sigma model result of as we will show below. 

Diagonal approximation and sigma model. For energy 
windows of size iVxh > VAk na , corresponding to r^r in 
Eq. . © smaller than the time of flight through the sys- 
tem, i.e. the ballistic Thouless time rxh, all contributions 
beyond the direct path are damped out. In this univer- 
sal regime the RWM predictions (consistent with RMT) 
are given by substitution of the direct path contribution 
to Ri,j (first line in Eq. JHJ) into Eqs. (@J or ©. It 
has been shown that the Gaussian and sigma-model re- 
sults coincide at this universal level 0, • For smaller 
energy windows, however, system-dependent deviations 
appear, i.e., T = F RWM + F SYS . In the present Gaus- 
sian approach deviations J^ G - SYS from universality are 
obtained by substitution of the non-direct contribution 
to Rij (second line of Eq. (JSJ) into Eq. (0. As we see, 
j:G,Sys cons i s t s of coherent single and double sums over 
non-direct classical paths of increasing length, while de- 
viations F a > SYS from universality in the BaM are ex- 
pressed through a purely classical object, the ballistic 
propagator |ll|. It is by no means clear whether the two 
approaches for !F SYS are consistent, and we address this 
fundamental issue now. 

The frequency-dependent ballistic propagator DV (w) 
is constructed by projecting the resolvent of the classi- 
cal Liouville equation, {H(r,p),p} — iwp (where {. . .} 
is the Poisson bracket), onto configuration space at en- 
ergy E na = p\ a j2m. The first step in unifying the lo- 
cal Gaussian and sigma-model approaches consists in ex- 
pressing IP'- 7 (w) through a sum over non-zero paths |2lj . 
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bution IIo from direct paths (set to zero when r\ 



contains £>° 4 
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The related energy-averaged, smoothed version reads 



ip j h = n; j w+Vfl 



7i,j 



TN 



(10) 



with IP J (u;) ~ rP J (u>) for tn ^> rxh- Consider now 
terms in the double sums in Eq. © where the classical 
paths involved join the same points, namely ff — rf = Ti 



The expression R % ^R^ 



is a double 
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sum over terms oscillatory in [S 7i } (En a ) ± {E^^ )]/h. 
When the energy difference hw a ,p — A(n^ — n a ) is clas- 
sically small, each trajectory in R % ^ is a smooth defor- 
mation of a corresponding trajectory in R)£ , and we can orbits where studied in |26 
expand S yi:j (En ) ~ S~ li j (E na ) + Tiw a .j3T~ li j . Collect- 
ing the non-oscillatory terms in the action difference of 
non-direct paths we obtain the so-called diagonal part 



T>i,j u l J 



SYS 
diag 



1 IT i- 1 "' J 



(11) 

Substitution of Eqs. (|10lll|) in Eq. (J5J shows that the 
diagonal part of p G > SYS is of order h e s and given by 



a linear combination of smoothed ballistic propagators 
(with universal coefficients Tf^ k l , J\f k ;) with the di- 
rect path contribution excluded. Using the general for- 
mulas (|(jl7fl the diagonal contribution to F G < SYS is read- 
ily calculated, giving in the limit tn ^> iTh exactly 
the various specific BctM results^^] available in the lit- 
erature (moments of the wavefunction, distribution of 
intensities, two-energies four-point correlations and 
two-point intensity distributions (20j). Hence we con- 
clude that to leading order in the deviation from uni- 
versality, the BaM corresponds to a Gaussian theory in 
diagonal approximation. For instance, Eq. © yields 
Zi/(27rTO 2 )[exp(-i/2)/V27rtK3-6/j + /; 2 ) for the prefactor 
of the ballistic propagator in the intensity distribution 
T G (t), in perfect agreement with supersymmetry [rH ]. 

Here several remarks are due: (i) Possibly most impor- 
tantly, our method provides also the general leading- 
order deviation from universality for the body of averages 
in disordered systems in the metallic regime by simply 
replacing IF (w) by the diffusive propagator; (ii) there 
is not a single chaotic system where the exact ballistic 
propagator IF^ (w) is known; hence the use of few classi- 
cal paths to construct the smoothed version (|10|l makes 
the Gaussian theory more accessible for practical calcu- 
lations; (iii) the direct path contribution, known to be 
counted twice as an artifact in the sigma model calcu- 
lations is correctly incorporated in the present ap- 
proach; (iv) our results can be easily generalized to the 
case of broken time-reversal symmetry (by taking each v a 
in Eq. as a complex vector with independent real and 
imaginary parts) and to the case of smooth potentials by 
using a sum over classical paths instead of Eq. JSJ. 

Statistics and oscillations of Coulomb Blockade trans- 
mission peaks. As a further application of our approach 
we consider transport through a quantum dot weakly 
coupled to two leads. In this Coulomb Blockade regime, 
characterized by the mean resonance width T <C k^T <C 
A, transport is mediated by resonant tunneling with cor- 
responding distinct conductance peaks. These Coulomb 
Blockade peaks and the fluctuations of their heights have 
been prominent objects of experimental studies [22|. The 
universal contribution to the conductance distribution, 
derived in RMT [2^ for the case of one-channel leads, 
was extended to the multichannel case for ballistic quan- 
tum dots using the RWM 24j and for disordered systems 
using the sigma model |2f | , while effects due to periodic 

However, for disordered dots 
in the Coulomb Blockade regime we are not aware of any 
prediction beyond universality. Moreover, non-universal 
oscillatory effects in ballistic dots were presented in , 
but the first-order theory used there fails to reproduce 
numerical results for the case of asymmetric leads |2gj . 

Here we apply our approach to these two problems. 
To this end we consider leads supporting one channel 
each, connected to the dot at positions fi (i = 1,2) 
with equal coupling strength a, generalizations to more 



4 



channels are straight forward 0|. Following [23I l24| . 
the height of the n Q -th conductance peak is given by 
G na = {e 2 /h)(aA)/(2Trk B T)g na where 



9n 



■kA 



V'L (n) + ^L(^)' 



(12) 



which fluctuates with n Q . Applying our general Eq. © 
we obtain the distribution of conductances as P(g) = 
(27r 5 )- 1 / 2 exp [-g/2] [l + 5PM(g) + 5P {2 \g)] with first 
and second order deviations from universality given by 



5P {1) {g) 
5P {2) {g) 



R 



2:2 



ql(g) (R^r + (R 2 ' 2 f 



(13) 



+ qf\g)R^R^ + q^(g)(R^) 



(2), 



Here the polynomial corrections to the RMT result, 
qW(x)(l/4)(x-l),q[ 2) (x) = (l/32)(a; 2 -9x+6), qf\x) = 

(l/16)(x 2 -3x), and qf ] = (1/8) (x 2 -7x + 4), are modu- 
lated by the single and double sums over classical paths. 

Using P(g), we now calculate the mean of the conduc- 
tance peak heights, (g) = 1 + <5_<7 ^ 1 ^ + 6g^ 2 \ First we 
discuss the leading term Sg^ = (1/2) (R 1 * 1 + R 2 ' 2 ) ~ 

1 /2 

0(h c ff ). In a billiard system with energy window Se ~ 
fo/T~Th we have R i,t ~ cos(kLi) where Li is the length of 
the shortest classical trajectory starting and ending at fi . 
The corresponding modulations of the conductance with 
frequencies Li,L 2 were already reported in [27I El ] . 

Our approach now enables us to go beyond the lead- 
ing semiclassical order and to provide the explicit ex- 
pression 5gW = -(3/16XR 1 ' 1 - i? 2 < 2 ) 2 - (l/4)(i?^ 2 ) 2 ~ 
0(h c s). As we see, Sg^ will be a combination of terms 
cos[k(Li ± Lj)]. We conclude that the second-order cal- 
culation is essential to understand the modulations with 
frequency L\ ± L2 numerically observed for asymmetric 
leads in Ref. [2^] . To our knowledge neither the first- nor 
the second-order oscillatory effect, present in Coulomb 
blockade physics, is accessible using the BerM. 

For diffusive dots in the metallic regime Sg^ is expo- 
nentially supressed and we just replace in the diagonal 
contribution to 5g( 2 ' the smoothed ballistic propagator 
IT 2J (u') by its diffusive counterpart H]iL( w ) 1° obtain 



Sgf 



dis 



3fi 
8irm'- 



Since n£ 



(0) > 0, we predict that the leading non- 
universal correction to the mean conductance peak height 
in weakly disordered quantum dots (with single-channel 
point contacts) is always negative and of order h e g. 

To summarize, we have used a local Gaussian theory 
for eigenfunction statistics to derive both smooth and 
oscillatory effects beyond (and including) the universal 
Random Wave Model. Smooth contributions are shown 
to give existing results of the Sigma Model, illuminating 



the connection between the two methods beyond univer- 
sality. In view of that, we use our approach to present 
new results for the conductance of diffusive quantum dots 
in the Coulomb Blockade regime. For the ballistic case, 
new oscillatory effects neglected both by the Sigma Model 
and previous semiclassical approaches are shown to de- 
scribe previously unexplained numerical results. 
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